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Abstract 

We manifest the origin of the wrong conclusion made by several groups of authors 
on the absence of Coulomb corrections to the cross section of the e + e~ pair production 
in ultrarelativistic heavy-ion collisions. The source of the mistake is connected with 
an incorrect passage to the limit in the expression for the cross section. When this 
error is eliminated, the Coulomb corrections do not vanish and agree with the results 
obtained within the Weizsacker- Williams approximation. 



RHIC and LHC projects initiated a set of recent publications on the e + e~ pair pro- 
duction in ultrarelativistic heavy-ion collisions. Using slightly different approaches, the 
authors of |], H, calculated the cross section of the process exactly in the parameters 
C(Za,b {Za,b being the charge numbers of the nuclei A and B, a is the fine-structure con- 
stant). In these papers the nuclei were treated as sources of the external field, and the 
amplitude was calculated at a fixed impact parameter of the nuclei. After that the cross 
section was obtained by the integration over the impact parameter. As a result, the conclu- 
sion was made that the exact cross section coincides with that calculated in the lowest order 
perturbation theory with respect to olZ>a,b (Born cross section). On the other hand, in the 
Weizsacker- Williams approximation with respect to one of the nuclei, the cross section of 
the process is proportional to the well-known cross section of the e + e~ pair production by 
a photon in a Coulomb field |3J] and, therefore, contains the Coulomb corrections. This 
obvious circumstance was observed in [[J, where the Coulomb corrections in the process 
under discussion were calculated. Though the existence of the Coulomb corrections is out 
of doubt, the source of the disagreement between the results was not revealed so far. This 
question is important from the theoretical point of view, since the approach developed in 
IHI, H, |j is used now in QCD. In the present paper we present the solution of this puzzle. 

Let the ultrarelativistic nuclei A and B move in the positive and negative directions 
of the z axis, respectively. Then the expression for the cross section of the e + e~ pair 
production, obtained in |], ^ ||, reads 

MM = U(p) «(k-Px)+7om 

—p+q- — (k — PjJ 2 — m 2 + xe 



_, , -ol (k - qx) + 7o™ , . 

'P-Q+ ~ (k — q_ij — mr + te 



Here p and e v (q and e q ) are the momentum and energy of the electron (positron), u(p) and 
u{—q) are positive- and negative-energy Dirac spinors, o: = 7°7, 7± = 7° ± , y z , 7 M are the 
Dirac matrices, p± = e p ± p z , q± = e q ± q z , m is the electron mass, k is a two-dimensional 
vector lying in the xy plane, and the function F(A.) is proportional to the electron eikonal 
scattering amplitude in the potential V(r) of the corresponding nucleus: 

oo 

F(A) = J d 2 pexp[-zpA] {expHx(p)] - 1} , X (p) = J dzV(z,p) . (2) 

— oo 

For the potential V(r) = V c (r) = —Za/r, the integral in x(p) becomes divergent and 
requires a regularization. This regularization can be made by using the potential V(r) = 
— Zaexp(— r/a)/r. Performing the integration in (0), and taking the limit a — > oo at fixed 
one obtains (up to the constant phase depending on a): 



T(l + iZa) VA 

Actually, to obtain this result one can use any regularization of the phase x(p) f° r which 
X{p) -> at p -»• oo. Since |^(A)| 2 = (4vrZa/A 2 ) 2 oc Z 2 , then the substitution (|) into 
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(H) would lead to the wrong conclusion [l], ||, |3| that the exact cross section coincides with 
the Born result. Let us show that, in order to obtain the Coulomb corrections in (|1|), it is 
necessary first to take the integral over k using the functions F(A) with the regularized 
phase and then remove the regularization. 

Consider the integral 

G = /(0 fc2 ( |F(k)|2 - |F ° (k)|2 ) ' (4) 

where F°(A) = —i J dpexp(—iA.p)x{p) is the first term of the expansion of F(A.) with 
respect to the potential. For F = T and, correspondingly, F° = JF° = AinZa/A 2 , the 
integrand in (§) vanishes. Let us show that the integral G is not equal to zero for the 
regularized F and is independent of the regularization method, if V(r) — ► —Za/r at r — ► 
( when x(p) ~ > 2Zaln(p) + const at p — > 0). For the sake of simplicity, we present the 
proof of this statement for a spherically symmetric potential V(r). Taking the integral in 
(§) over the angle of p, and integrating by parts over p, we obtain the following expression 
for F: 

oo 

F(k) = — J dppJ 1 (kp)x'{p)exp(-ix{p)) , (5) 
o 

where J±(x) is the Bessel function. The function F°(k) can be obtained from fl5|) by omitting 
the exponent in the integrand. Substituting (|5]) into (|]), and integrating over the angle of 
k, we find 

oo oo oo 

G = 2n J dqq J J dpi dp 2 PiP2Ji(kp 1 )J 1 (kp 2 ) x 



x x'(pi)x'(p2) {exp[-zx(pi) + ix(p2)] - 1} • (6) 

If one changes naively the order of integration over k and and takes the integral over 
k, using the relation 

oo 

/ dk kJ 1 (kp 1 )J 1 (kp 2 ) = 5(pi - p 2 ) , 

then, after the integration over pi, the result will be zero. To demonstrate that the change 
of the integration order in (^) is invalid, we restrict the upper limit of the integral over k by 
Q. After that one can change the order of integration in triple integral in (f|). Integrating 
over k, we obtain 

oo oo 

G = 2tt dp 1 dp 2 ^-^[p 2 J (Qp 2 )J 1 (Qpi) - piJ (Qpi)Ji(Qp2)} x 

ri 

xx'(Pi)x'(P2) {exp[-zx(pi) + ix{p2)] ~ 1} ■ (7) 

Substituting p\ y2 — > pi t2 /Q, and taking the limit Q — > oo with the use of the asymptotics 
of x, we find 



oo oo ^ 

G = 8n(Za) 2 I I dp l dp 2 — 2 [p2Jo{p2)Ji{pi) ~ Pi J o(pi) J i{p2)} * 

J J jOi — Oo 




Pi ~ P2 

I } ■ (8) 
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Making the change of variables — Texp(±t/4), and integrating over r, we obtain the 
non-zero result for the quantity G: 



DO 



,, / , cos( Zat) — 1 

G = 8tt (Za 2 / dt — 9 

' exp(t) - 1 y 1 



= -87r(Za)'[Re^(l + iZa) + C) = -^{Zaf f{Za) , 

where C is the Euler constant, ip(x) = d\nT(x)/dx. Thus, we come to the remarkable 
statement: although the main contribution to the integral in comes from the region of 
small k, where |-F(k)| differs from |.F(k)| = 4irZa/k 2 and depends on the regularization 
parameters (the radius of screening), nevertheless, the integral G itself is a universal func- 
tion of Za. Note that the integral (H) appears in the theory of multiple scattering (see || 
where the approximate formula for this integral was obtained). 

Now it is clear, how to derive the Coulomb corrections starting from the expression 
(H). Let us calculate the Coulomb corrections related to the nucleus B (the contribution 
of the higher order perturbation theory with respect to the parameter Zbol)- For this 
purpose one should replace in flU) the functions |-Fb| 2 and \Fa\ 2 with |-Fb| 2 — \F^\ 2 and 
|JF^| 2 , respectively, keeping the regularization in the functions Fb and F B . The main 
contribution to the integral is given by the region of small k. Therefore, we can neglect k 
in the argument of T\ and expand the matrix element M. at small k: 

MW » kL , L = SM ( g (->-/»+- 2±M + ^Jl± - } „(_,) . (10 ) 

I (p_ + q_) (p_ + q_y J 



Using (P) and (|10|) , and performing the summation over electron and positron polarizations, 
we obtain the following expression for the Coulomb corrections related to the nucleus B\ 

c 2G B d 3 pd 3 q |^(p ± + q±)| 2 / , , 2 2(p ± g + g_ + q ± p + p_) 2 | 

Here denotes the function G in at Z — Zb- The Coulomb corrections related to the 
nucleus A can be obtained from (|TT| ) by the substitution and the replacement of 

indices — ^> +. 

It is necessary to note the following circumstance. Actually, in the expansion over Z^a 
and Z B oi of the differential cross section da/dpdq in ([[]) , only the lowest (Born) term is 
correct. As for the higher order terms in (|l|) (Coulomb corrections), they give the correct 
result only after the integration over the directions of the positron (electron) momentum. 
This is due to the fact that the asymptotic form of the wave functions in |], |j, [| corresponds 
to the problem of scattering, but not to the problem of pair production. If one calculates the 
cross section integrated over the direction of q, then, due to the completeness relation, it is 
possible to replace the set of functions containing in asymptotics the converging spherical 



wave with the set of functions containing the diverging spherical wave. Thus, (|TT|) should 
be integrated over the angles of q (p). The same trick was made at the recalculation of the 
bremsstrahlung cross section integrated over the photon momentum from the cross section 
of pair photoproduction integrated over the positron momentum 0. It explains why the 
Coulomb corrections ( fTT| ) are given by the region of small k, while at the calculation of 
the Coulomb correction using the wave functions with the correct asymptotic behavior the 
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main contribution would come from the region k ~ m. The same situation occurs at the 
calculation of bremsstrahlung and pair photoproduction cross sections, where the Coulomb 
corrections come from different regions of momentum transfers. 

Let us calculate within the logarithmic accuracy the Coulomb corrections to the cross 
section da/de p de q at e PA ^> m. At the integration over the transverse momenta the main 
contribution comes from the region A = |p_i_ + q±| <C pi, ?i ~ m. The integral over 
A requires regularization at A — > 0. It is obvious that the lower limit of integration 
over A coincides with that in the Weizsacker- Williams method. In the rest frame of the 
nucleus B it has the form A min = (e° + e°)/7, where e^ q are the energies of the electron 
and positron, 7 is the Lorentz factor of the nucleus A in this frame. In the laboratory 
frame, where the nuclei A and B have the Lorentz factors y A and jb, respectively, one has 
A min = (p + + q + )/j A . Using this cutoff, we obtain 



2 2 

m m 
In —75 h m 



2min 



A 2 1 A 2 

(12) 

The sum of logarithms in this formula corresponds to the contributions of two kinematic 
regions: p z , q z > and p z , q z < 0. In the first case Ai min = {e p + e q )/y A , and the 
corresponding term in (|12[) is valid at m < e Pt(] <C my A . In the second case A 2min = 
m 2 /(e p + e q )yA, and the corresponding term is valid at m C e PtQ m^B- Performing the 
integration over e Ptq in the regions indicated, one has 

28 

0% = -^(Z A a) 2 (Z B a) 2 f(Z B a) \n 2 ( lAlB ) ■ (13) 

The formulas ( |12"D and (p~3f ) can be easily obtained in the Weizsacker- Williams approxima- 
tion using the well-known result for the exact in Za pair photoproduction cross section in 
the field of a nucleus. They coincide with the result of (see also 0). It also follows from 
the Weizsacker- Williams method that the contribution of the terms, containing the higher 
orders of Z A and Zb simultaneously, can be neglected within our accuracy. 

In papers |T], 0, [| the amplitude of e + e~ pair production was obtained at fixed impact 
parameter between the nuclei. Using this amplitude, it is possible to represent the Coulomb 
corrections related to the nucleus B as the integral over the impact parameter: 

^ - T^S/ d2 "//SS exp[i(ki - k2)plAi(ki)Ai * (k2)x (i4) 



X 



>B(k x )^(ka) - ^(kO^Cka)! ^(qx + P± ~ kO.F°*(q ± + p ± - k 2 



Again, changing the order of integration would lead to zero result. Indeed, taking the 
integral over p first, we get the factor S(ki — k 2 ) in the integrand, and, therefore, the 
integral over ki vanishes due to the relation {J-'bI 2 = I-^bI 2 - Let us demonstrate that, 
similar to the case of the integral (|j) calculation, the change of the integration order in 
(|H|) is incorrect, and the result ([□]) also follows from (|Tj). For this purpose, we restrict 
the region of integration over p by the condition p < R. After that it is possible to change 
the order of integration and take the integral over p. Then the main contribution to the 
integral over k 12 comes from the region /c 12 < 1/R. Since we are going to take the limit 
R — > 00, we can replace .M(ki j2 ) with k ij2 L and neglect ki j2 in J r A (q± + Pj_ — ki j2 ). Then, 
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we have 



^Z7TJ tp£q 



OCSX) 








x [k 2 RJ (k 2 R)Ji(k 1 R) - k 1 RJ (k 1 R)J 1 (k 2 R)} . 

Comparing the expression for the function Gb with we see that Gb = Gb- After the 
summation over the electron and positron polarizations the formula ( f[5|) comes into (11). 
Note that the expression ( |TTj) can be obtained directly from (|I4]) by taking the integral over 
k 12 in the region /c 12 < |pj_ + q±| and then integrating over p in the infinite limits. 

If the impact parameter p is restricted by the beam transverse size R®, then it follows 
from the above consideration that the effect of the finite size appears when we can not 
neglect k ~ 1/Rq in the argument of T\ in comparison with | q± + Pj_ | . This is equivalent 
to the condition R -C 1/A min ~ ■JaIb/^i- 

Thus, the method developed in |I|, |2|, |3| can be used for the calculation of the Coulomb 
corrections to the e + e~ pair production cross section integrated over the direction of the 
positron (electron) momentum. Its careful application leads to the correct result. 

We are grateful to V.M. Katkov and V.M. Strakhovenko for useful discussions. 
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